Comparative and phylogenetic analyses of Loranthaceae plastomes provide insights into the evolutionary trajectories of plastome degradation in hemiparasitic plants

Background The lifestyle transition from autotrophy to heterotrophy often leads to extensive degradation of plastomes in parasitic plants, while the evolutionary trajectories of plastome degradation associated with parasitism in hemiparasitic plants remain poorly understood. In this study, phylogeny-oriented comparative analyses were conducted to investigate whether obligate Loranthaceae stem-parasites experienced higher degrees of plastome degradation than closely related facultative root-parasites and to explore the potential evolutionary events that triggered the ‘domino effect’ in plastome degradation of hemiparasitic plants. Results Through phylogeny-oriented comparative analyses, the results indicate that Loranthaceae hemiparasites have undergone varying degrees of plastome degradation as they evolved towards a heterotrophic lifestyle. Compared to closely related facultative root-parasites, all obligate stem-parasites exhibited an elevated degree plastome degradation, characterized by increased downsizing, gene loss, and pseudogenization, thereby providing empirical evidence supporting the theoretical expectation that evolution from facultative parasitism to obligate parasitism may result in a higher degree of plastome degradation in hemiparasites. Along with infra-familial divergence in Loranthaceae, several lineage-specific gene loss/pseudogenization events occurred at deep nodes, whereas further independent gene loss/pseudogenization events were observed in shallow branches. Conclusions The findings suggest that in addition to the increasing levels of nutritional reliance on host plants, cladogenesis can be considered as another pivotal evolutionary event triggering the ‘domino effect’ in plastome degradation of hemiparasitic plants. These findings provide new insights into the evolutionary trajectory of plastome degradation in hemiparasitic plants. Supplementary Information The online version contains supplementary material available at 10.1186/s12870-024-05094-5.


Background
Parasitic plants rely on their host to alleviated competition with other plant species for essential resources (e.g., light, water, soil nutrients), thereby enhancing their adaptive capacity in unfavorable environmental conditions; consequently, parasitism is generally perceived as a life strategy that evolves under natural selection and environmental pressures [1].To date, approximately 4,750 parasitic species have been documented in angiosperms and it is estimated that the lifestyle transition from photosynthetic autotrophy to parasitism may have independently evolved at least 12 or 13 times in the angiosperm tree of life [1,2].Based on their photosynthetic capacity, parasitic plants are divided into two categories: hemiparasites and holoparasites [3][4][5].Hemiparasites retain varying degrees of photosynthetic capacity while obtaining nutrients (primarily water and mineral elements) from their hosts.Contrarily, holoparasites exhaustively lose their photosynthetic capacity, with severely or entirely degraded leaves; therefore, they must completely rely on the host plants for energy and nutrients [3,5,6].
Chloroplasts, which originated from the endosymbiosis of photosynthetic autotrophic cyanobacteria within primitive eukaryotic cells, serve as vital centers for energy acquisition and carbohydrate production in plants [7,8].Due to their distinct origin, chloroplasts possess autonomous genomes and transcription/translation systems that synergistically interact with the nuclear genomes for plastid protein synthesis [8,9].Although the chloroplast genome (plastome) is highly conserved in most photosynthetic autotrophic angiosperms in terms of genome size, structure, and gene content, the lifestyle transition from autotrophy to heterotrophy might relax the purifying selection pressures that play a crucial role in maintaining the stability of plastomes [10], resulting in varying degrees of gene loss, pseudogenization, size reduction and structural rearrangement in parasitic plastomes [10][11][12][13][14][15][16].Previous studies have demonstrated that different levels of reliance on host plants for energy and nutrient requirement can exert various intensities of selective pressures on the plastomes of parasitic plants; as a result, holoparasites exhibit higher levels of plastome degradation compared to their hemiparasitic relatives [13][14][15][16][17][18][19].
Hemiparasitic plants are further classified into facultative root-parasites and obligate stem-parasites based on their feeding modes and nutritional dependence on the hosts [20].Considering the relatively higher degrees of trophic reliance of obligate stem-parasites on their hosts [21], purifying selection of plastid genes tends to be gradually relaxed during the transition from facultative parasitism to obligate parasitism.Therefore, obligate stem-hemiparasites are expected to have higher levels of plastome degradation than facultative root-hemiparasites [14,22].Surprisingly, this theoretical prediction has not yet been empirically validated, as previous comparative studies on the plastomes of Santalales hemiparasites revealed that obligate stem-parasites did not consistently exhibit a higher level of gene loss or pseudogenization than facultative root-parasites [23-25,].Additionally, a previous study suggested that plastome degradation in hemiparasitic plants may follow a 'domino effect' , where initial losses of plastid genes trigger chain reactions leading to subsequent gene losses [23].However, the evolutionary mechanisms underlying the enigmatic 'domino effect' of plastome degradation in hemiparasitic plants remain elusive.These unresolved issues leave critical gaps in the exploration of evolutionary trajectories of plastome degradation in hemiparasitic plants.
Previous studies have demonstrated multiple independent evolutionary transitions from facultative rootparasitism to obligate stem-parasitism in Santalales [1,24,25], and proposed that plastome degradation in Santalales hemiparasites may have evolved in a lineagespecific manner [24][25][26].This suggests that the degree of plastome degradation may exhibit substantial variation among phylogenetically distant facultative root-parasites, potentially surpassing the differences in plastome degradation observed between closely related root-parasites and obligate stem-parasites.Consequently, previous studies [24,25,27] on the retrogressive evolution of hemiparasitic Santalales plastomes may be flawed, as comparative analyses of plastomes based solely on lifestyle variations (facultative root-parasitism vs. obligate stem-parasitism) without considering phylogenetic relationships could lead to biased inferences.
The hemiparasitic lineage Loranthaceae, which includes approximately 76 genera and more than 1,000 species, is the largest family within Santalales [20,28,29].Within Loranthaceae, three monotypic genera, Nuytsia, Atkinsonia, and Gaiadendron, are facultative root-parasites, while the remaining genera are obligate stem-parasites [20].Given the monophyletic origin of obligate stem-parasitism from facultative root-parasitism in Loranthaceae [20,29,30], this family provides an ideal system to investigate whether obligate stemparasites have a higher degree of plastome degradation than closely related facultative root-parasites.To gain a better understanding of the evolutionary trajectories of plastome degradation in hemiparastic plants, this study aims to (1) verify the theoretical prediction that the lifestyle transition from facultative root-parasitism to obligate stem-parasitism may lead to further degradations of hemiparasitic plastomes [14,27] and (2) explore the plausible events that could have triggered the 'domino effect' in the plastome degradation among Loranthaceae hemiparasites.To achieve these objectives, we sequenced the plastomes of 22 Loranthaceae hemiparasites and combined them with publicly available data, yielding an extensive dataset comprising the plastomes of 48 species from 14 genera within this family.Under the framework of plastome phylogeny, we analyzed plastome degradation in each clade to determine whether there was lineage-specific gene loss, pseudogenization, or plastome structural variation within Loranthaceae.Additionally, we compared the plastome features of obligate stem-parasites with those of facultative root-parasites to explore differences in the degree of plastome degradation.

Plant materials, illumina sequencing, plastome assembly and annotation
In total, 22 species from eight genera of Loranthaceae were collected from the wild (voucher specimens were deposited in the herbarium of Yunnan University and taxonomically identified by Dr. Min Deng), within voucher information being presented in Table S1.The genomic DNA of each sample was extracted from approximately 100 mg of silica gel-dried leaf tissues using a modified CTAB method [31].The purified genomic DNA was use to prepare a paired-end Illumina sequencing library (with an average insert size of approximately 350 bp) for each sample, employing the TruSeq DNA Sample Prep Kit (Illumina, Inc., San Diego, CA, USA) following the manufacturer's guidelines.Low-coverage genome sequencing was performed on an Illumina HiSeq 2500 platform to generate approximately four Gbp raw reads for each sample.
Raw Illumina reads were filtered using an NGS QC tool kit [32] to remove adaptors and reads with ambiguous bases.The complete plastome of each sample was de novo assembled using NOVOPlasty v.2.7.0 [33] with the k-mer size set at 31, and the large subunit of the RuBisCO gene (rbcL) of Taxillus sutchuenensis (MG999457) was utilized as the seed for plastome assembly.These assembled plastomes were annotated using PGA [34], and the boundaries of start/stop codons and introns/extrons of protein-coding genes were checked using Geneious v10.2 [35].The tRNAscan-SE 1.21 [36].was used to verify the transfer RNA (tRNA) genes using default parameters.All annotated plastomes were deposited in the NCBI GenBank database and their accession numbers are listed in Table S1.In addition to these newly sequenced plastomes, publicly available complete plastomes of 26 Loranthaceae species (Table S2) were obtained from the NCBI GenBank database for phylogenetic and comparative analyses.These publicly available plastomes were reannotated using the same method.

Phylogenetic analyses
In total, 48 species of Loranthaceae hemiparasites were sampled for phylogenetic analysis.Erythropalum scandens (Erythropalaceae), an autotrophic relative of Loranthaceae [1,28,30] was selected as the outgroup.Based on the plastome dataset, 55 plastid protein-coding genes (PCGs) commonly shared by these species were extracted from each plastome using Geneious v10.2 [35].These PCGs were aligned separately using MAFFT v7.3.1 [37], and each CDS matrix was concatenated into a supermatrix (3,8153 bp) using Phylosuite [38].Based on the concatenated matrix, maximum likelihood (ML) and Bayesian inference (BI) approaches were employed to infer the intrafamilial phylogenetic relationships of Loranthaceae.ML analysis was implemented by IQ-TREE V2.1.2[39,40] using UFBoot2 with the best-fit model TVM + F + R2 determined by Modelfinder [41]; 1,000 replications were adopted to calculate the standard bootstrap percentage for each node.MrBayes v3.2.7a [42] was used for the BI analysis.According to the Akaike information criterion (AIC), the best-fitting partitioning schemes and models were selected for each gene with default values using PartitionFinder v2.1.1 [43]. .The Markov chain Monte Carlo algorithm was run for five million generations with every 5,000 generations for tree sampling.Trees resulting from the first 25% of generations were discarded as "burn-in", and the remaining trees were used to build the 50% majority-rule consensus tree and estimate the posterior probability values (PP).FigTree 1.4.6 [44] was used to present and edit the phylogenetic trees.

Comparative analyses of Loranthaceae plastomes
Genomic features, including genome size, gene content, and GC percentage of the whole plastome, LSC regions, SSC regions, and IR regions, were compared according to the sequence and annotation using Geneious v10.2 [35].The boundaries of LSC, SSC, and IRs of each plastome were compared to check IR expansion/contraction in Loranthaceae plastomes using Geneious v10.2 [35] by comparison with their autotrophic relative, Erythropalum scandens.The microsyntenic structure of Loranthaceae plastomes was investigated using the multiple genome alignment tool, Mauve v2.3.1 [45], after removing one copy of the IR region for each plastome.Gene loss and pseudogenization events detected in each plastome were mapped to a phylogenetic tree to trace the evolutionary trajectory of plastome degradation associated with the evolution of the hemiparasitic lifestyle in Loranthaceae.

Phylogeny-oriented comparative analysis of Loranthaceae plastomes
Except for a slight fragment translocation with a sequence length of 25 bp in Phyllodemis delavayi, microsynthetic analysis showed that no additional structural rearrangements occurred in these Loranthaceae plastomes (Fig. 2).Despite obvious variations in genome size (Table 1), these Loranthaceae plastomes exhibited a high level of conservatism with respect to their LSC/ IRa, IRa/SSC, SSC/IRb, and IRb/LSC junctions (Fig. 3).Although Loranthaceae hemiparasites shared same SSC/ IRb and IRb/LSC boundaries with the autotrophic relative Erythropalum sandens, the LSC/IRa boundaries of Loranthaceae plastomes moved into the rpl2 gene (Fig. 3).Within Loranthaceae, the SSC/IRa boundary of the facultative root-parasite Nuytsia floribunda moved into the rpl32 gene, whereas the SSC/IRa boundaries of all obligate stem-parasites were further contracted into the trnL-UAG gene (Fig. 3).
To trace the evolutionary trajectory of plastome degeneration associated with the formation of hemiparasitic lifestyle in Loranthaceae, the gene loss and pseudogenization events detected in each species were mapped to a phylogenetic tree (Fig. 4).As ycf15 was annotated as a pseudogene in all analyzed plastomes, this mutation may have occurred in the most recent common ancestor (MRCA) of Loranthaceae and Erythropalaceae (Fig. 4).Contrary to their autotrophic relative, Erythropalum scandens, all plastid NADH dehydrogenase-like complex genes, except for ndhB, two ribosomal protein genes (rps15 and rps16), and two tRNA genes (trnG-UCC and trnV-UAC), were absent from Loranthaceae plastomes, which may represent a molecular synapomorphy of the family (Fig. 4).Compared to the facultative root-parasite Nuytsia floribunda, additional deletion of rpl32 and ndhB genes might have taken place in the stem lineage ancestor of obligate stem-parasites (Fig. 4).Among Loranthaceae obligate stem-parasites, further loss of the infA gene is commonly observed in the tribe Elytranthaea (Fig. 4).Within the tribe Loranthaea, the pseudogenization of pasI and infA might have occurred in the stem lineage ancestor of the subtribe Loranthinae, and the losses of trnA-UGC and trnL-GAU likely occurred in the MRCA of the subtribes Dendrophthoinae (except for Helixanthera parasitica), Dendrophthoinae, Amyeminae, Tapinanthinae, and Emelianthinae (Fig. 4).Furthermore, the pseudogenization of plastid rpl16 gene could have taken place in the MRCA of the subtribes Scurrulinae and Helixanthera parasitica (Fig. 4).

Characterization of plastome degradation in Loranthaceae hemiparasites
A growing body of evidence suggests that the lifestyle transition from autotrophy to heterotrophy often results in rampant plastome degradation in parasitic plants, including the physical and functional loss of plastid genes, genome downsizing, and structural mutations [13,14,16,23].Consistent with this, the current study found that the plastomes of 48 Loranthaceae hemiparasites exhibited varying degrees of plastome size shrinkage, structural alteration, gene loss, and pseudogenization compared to their autotrophic relative, Erythropalum scandens.The findings indicate that although the hemiparasitic plants enable photosynthesis, their plastomes have experienced diverse modifications owing to the evolution of heterotrophic lifestyles.
Compared to Erythropalum scandens, gene loss and pseudogenization events occurring in the plastomes of Loranthaceae hemiparasites involved the following PCGs: ndhA-K, infA, rpl16, rpl32, rpl33, rpl36, rps15, rps16, psaI, psbZ, matK, ycf1, ycf2, ycf15, and rpoC2, as well as six tRNAs.Among them, only three PCGs (rps16, ycf2, and ycf15) and four tRNAs (trnA-UCC, trnH-GUG, trnI-GAD, and trnV-GAC) were located in the IR regions, whereas the remaining genes were located in the LSC and SSC regions.Although the plastomes of Loranthiaceae hemiparasites possessed lower GC content in the LSC and SSC regions than Erythropalum scandens, there were slight differences in the GC content in the IR regions.The higher levels of gene loss/pseudogenization and lower GC content in the LSC and SSC regions, in contrast to the lower levels of gene loss/pseudogenization and higher GC content in the IR regions observed in the plastomes of Loranthaceae hemiparasites, provide robust evidence for the theoretical prediction that loss/pseudogenization of plastid genes in parasitic plants may occur simultaneously with a decrease in GC content [16].
Previous studies have demonstrated that parasitic plants with low GC content in their plastomes tend to accumulate inversions and structural mutations around the IR regions [13,14].Although varying degrees of plastid gene loss occurred in Loranthiaceae hemiparasites, significant synteny was observed among their plastomes, suggesting that they did not undergo dynamic structural rearrangement after the formation of hemiparasitic lifestyle.Notably, compared with their autotropic relative, Erythropalum scandens (38.0%), the GC content of Loranthaceae plastome (35.9-37.8%)was only slightly decreased.The findings provide additional empirical evidence that a dramatic decrease in GC content is more Fig. 4 Comparison of gene content among Lorantahceae plastomes and Erythropalum scamdens plastome.Genes mentioned above branches indicate loss of genes, whereas under branches indicate pseudonization of genes.Red squares: intact genes; blue squares: deleted genes; yellow squares: pseudogenes likely to trigger the structural rearrangement of the parasitic plastomes [13,14].Due to a higher degree of plastid gene loss and reduction in GC content observed in the plastomes of holoplastic plants compared to hemiparasitic plants, the plastomes of holoparasites are more prone to independent structural rearrangement, as previous study demonstrated [46][47][48][49].
The gene contents of Loranthaceae plastomes and their autotrophic relative were compared based on the phylogenetic tree topologies.The results indicated that a total of 10 ndh genes (ndhA, ndhC-K), two ribosomal protein genes (rps15 and rps16), and two tRNA genes (trnG-UCC and trnV-UAC) might have been deleted in the stem lineage ancestor of Loranthaceae.The deleted ndh genes encode 10 subunits of the plastid NDH complex that mediates photosystem I cyclic electron transport and facilitates chlororespiration in plant cells [50].The functional and physical loss of most plastid ndh genes has been commonly observed in a wide spectrum of parasitic angiosperm lineages and is regarded as an early response to the evolution of heterotrophic lifestyle [13,14,16].Notably, this mutation has also been detected in a wide spectrum of photoautotrophic angiosperms [48,[51][52][53][54], and Frailey et al. proposed that the plastid ndh genes may have undergone negative selection in these photoautotrophic lineages [49].Additionally, Lin et al. observed that several ndh genes encoded by the chloroplast and nuclear genomes may have been lost concomitantly in some epiphytic autotrophic orchids and proposed that this mutation may increase the possibility of evolving a heterotrophic lifestyle [19].On this basis, we speculate that the massive reduction in the NDH pathway in the stem lineage ancestor of Loranthaceae might have played a pivotal role in triggering the transition from photoautotrophic to hemiparasitic lifestyles.
In addition to the loss of 10 plastid ndh loci, the deletion of two plastid ribosomal protein-encoding genes (rps15 and rps16) was inferred to have occurred in the stem lineage ancestor of Loranthaceae, and further loss/ pseudogenization of other plastid ribosomal proteinencoding genes (rpl16, rpl32, rpl33, and rpl36) was commonly observed in Loranthaceae obligate stem-parasites.Notably, loss/pseudogenization of these plastid ribosomal protein-encoding genes has been detected not only in parasitic plants [55,56], but also in a wide range of autotrophic angiosperms [57][58][59][60].Although the rps15, rpl33, and rpl36 genes are not essential for chloroplast gene translation, the remaining genes play essential roles [61].Hence, the loss of these plastid ribosomal proteinencoding genes in Loranthaceae may be compensated for by other plastid rpl/rps genes or by nuclear-encoded rpl/rps genes.
Plastid infA is another commonly reduced gene in the plastomes of Loranthaceae hemiparasites.While this gene was intact and retained in the facultative root-parasite Moquiniella rubra, the remaining Loranthaceae species (obligate stem-parasites) exhibited either loss or pseudogenization of this gene.This suggests that, after the evolution of a hemiparasitic lifestyle, plastid infA gene was recurrently deleted or pseudogenized in Loranthaceae obligate stem-parasites.Similarly, parallel losses of plastid infA genes have been observed in Santalales [24][25][26]62] and a majority of holoparasitic plants [14,16], as well as in diverse lineages of photoautotrophic angiosperms [16,63].Previous studies proposed that the plastid infA gene may have been constantly transferred to the nucleus [64].Accordingly, horizontal transfer of the infA gene from plastomes to nuclear genomes may have occurred independently in Loranthaceae obligate stem-parasites.
In addition to the abovementioned protein-encoding genes, the deletion of six plastid tRNA genes was observed in the plastomes of Loranthaceae hemiparasites.It is well known that some heterotrophic plants generally contain only a fraction of tRNAs [14,56,65].Despite the essential roles of tRNAs in plastid translation, the partial loss of these genes in plants may not represent a lethal mutation.Because of their relatively small size, tRNAs are easily transferred from the cytosol to plastid organelles [16], and photosynthetic plants most likely possess a specific transport mechanism that transports tRNAs to plastids from the cytosol [60,65].Accordingly, the loss of these tRNAs in Loranthaceae hemiparasites may not have a serious impact on their survival and is therefore tolerable.

New insights into evolutionary trajectories of plastome degradation in hemiparasitic plants
Along with the lifestyle transition from autotrophy to heterotrophy, parasitic plants developed numerous survival strategies to adapt to various environments, which might have led to diverse plastome degradation trajectories [15-17, 56, 66].Among hemiparasitic plants, facultative root-parasites show a slight degree of trophic reliance on the host compared to obligate stem-parasites, which require host support water and mineral elements to complete their life cycle after germination [20,21].Theoretically, the varying nutritional dependence of facultative root-parasites and obligate stem-parasites on their hosts may influence the reductive evolution of their plastomes [14,22].Although this theoretical prediction seems plausible, it has not been validated by empirical studies, leaving a critical gap in the exploration of the evolutionary trajectory of plastome degradation in hemiparasitic plants.
In the current study, the plastome phylogeny resolved all obligate stem-parasties as a well-supported monophyletic lineage sister to the facultative root-parasite, Nuytia floribunda, confirming the monophyletic origin of obligate stem-parasitism from facultative root-parasitism in Loranthaceae proposed by previous studies [20,29,30].Compared with their autotrophic relative, Erythropalum scandens, the plastomes of Loranthaceae hemiparasites exhibit various degrees of plastome downsizing and gene loss/pseudogenization.Among them, the plastome size of the early diverged facultative root-parasite, Nuytsia floribunda, was 139,027 bp, retaining 98 functional genes, which represents the largest plastome size and highest gene content in the Loranthaceae hemiparasites sampled in this study.Contrarily, the plastomes of all obligate stem parasites possessed relatively small genome sizes (115,635-129,570 bp) and retained relatively few functional plastid genes.Comparative analyses of Loranthaceae plastomes under the phylogenetic framework also showed that the stem lineage ancestor of obligate stem-parasites might have undergone more gene deletions (losses of rpl32 and ndhB genes) than the facultative root-parasite Nuytsia floribunda (pseudogenization of ndhB gene).Taken together, these results provide empirical evidence supporting the theoretical expectation that the evolution of obligate stem-parasitism from facultative root-parasitism might have caused higher levels of plastome degradation in hemiparasites [14,22].
Nevertheless, the taxonomic sampling in this study is still subject to certain limitations.Specifically, the Loranthaceae family has been documented to include three monotypic genera (Nuytsia, Atkinsonia, and Gaiadendron), each hosting a facultative root-parasite [20].Among them, Gaiadendron exhibits a closer phylogenetic relationship with the clade of obligate stem-parasites compared to Nuytsia and Atkinsonia [29].Given that this study only included one facultative root-parasite (Nuytia floribunda), comparing its plastome with those of obligate stem-parasites may not provide an accurate reflection of the evolutionary trajectory of plastome degradation associated with the lifestyle transition in hemiparasites.From this perspective, our data may not provide robust support for the theoretical expectation that the transition from facultative root-parasitism to obligate stem-parasitism could lead to increased levels of plastome degradation in hemiparasites [14,22].Therefore, the further validation is necessary to determine whether the lifestyle transitions resulting in increased nutritional reliance on host plants constitute pivotal evolutionary events that trigger the 'domino effect' in plastome degradation of hemiparasitic plants.
Within the clade of obligate stem-parasites, only one gene (infA), was lost in the plastomes of the early divergent tribe Elytranthaea.Within the tribe Loranthaea, further lineage-specific gene loss and pseudogenization occurred in the stem lineage ancestors of the subtribe Loranthinae, the MRCA of the subtribes Dendrophthoinae (except for Helixanthera parasitica), Dendrophthoinae, Amyeminae, Tapinanthinae, and Emelianthinae, as well as the MRCA of the tribes Scurrulinae and Helixanthera parasitica.Additionally, independent gene loss and pseudogenization were observed in some shallow nodes of the Loranthaea tribe.As a result, taxa within shallow branches exhibit a heightened degree of physical and functional loss of plastid genes compared to taxa within deep clades (e.g., the tribe Elytranthaea).This suggests that plastome degradation in Loranthaceae hemiparasites is likely a gradual process that intensifies with lineage and species divergence.Considering the occurrence of lineage-specific gene loss/pseudogenization at deep nodes and additional independent genes loss/pseudogenization in shallow branches, it is plausible to hypothesize that cladogenesis may act as an essential evolutionary force triggering the cascading degradation (the 'domino effect') of plastome degradation in hemiparasitic plants.

Phylogenetic implications
The infra-familial relationships of Loranthaceae recovered in this study were largely congruent with previous studies [20,29,30,[67][68][69] but received full support (BS = 100; PP = 1.00) for most nodes.Notably, the plastome phylogeny failed to resolve the genus Helixanthera as monophyletic, given that Helixanthera parasitica formed a sister relationship with the subtribe Scurrulinae clade, which is consistent with the findings of Su et al. [70] and Liu et al. [29].This relationship is supported by morphological evidence that the flowers of Helixanthera parasitica and the two genera (Scurrula and Taxillus) of the subtribe Scurrulinae are 5-merous, which is different from the 4-merous flowers of Helixanthera terrestris and Helixanthera sampsonii [29].Collectively, this line of evidence suggests that the current taxonomic delineation of the genus Helixanthera is likely problematic.To reasonably resolve this taxonomic issue, further studies based on multidisciplinary methods are required.
Consistent with a previous study [29], the plastome phylogeny showed that the two genera (Elytranthe and Macrosolen) of the tribe Elytranthaea were not monophyletic, as the species of the two genera were interwoven in the tree topology.Notably, both genera share similar vegetation, flower, and fruit morphology, and the diagnostic characteristics used to discriminate between them is that the inflorescences of Macrosolen have relatively larger bracts.Accordingly, Barlow argued that this trait may not be robust enough to distinguish these two genera from each other and proposed treating Macrosolen and Elytranthe as congeneric [71].The current study provides a phylogenomic evidence supporting this taxonomic proposal.

Conclusion
The plastome features of 48 Loranthaceae hemiparasites were characterized in this study.Under a phylogenetic framework, comparative analyses of these Loranthaceae plastomes revealed that all obligate stem-parasites exhibited a higher degree of plastid degeneration than the closely related facultative root-parasite, Nuytia floribunda, providing empirical evidence for the theoretical expectation that the evolution of obligate stem-parasitism from facultative root-parasitism may have caused a higher level of plastome degradation in hemiparasites [14,22].Phylogeny-oriented comparative analyses of Loranthaceae plastomes have also revealed that plastome degradation in Loranthaceae is likely a gradual process that intensifies with lineage and species divergence.Therefore, cladogenesis can be a key evolutionary force that triggered the 'domino effect' in plastome degradation of hemiparasitic plants.These findings provide new insights into the evolutionary trajectory of plastome degradation in hemiparasitic plants.

Table 1
Comparison of size and GC content (GC%) of complete plastomes, large single copy (LSC), small single copy (SSC), and inverted repeat (IR) regions, and gene content of Loranthaceae species with the autotrophic plant Erythropalum scandens